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^ ■ Abstract: We discuss the fundamental role of entanglement as the essential nonclas- 

. sical feature providing the computational speedup in the known quantum algorithms. 

O I We review the construction of the Fourier transform on an Abelian group and the 

Q I principles underlying the fast Fourier transform algorithm. We describe the imple- 

OO • mentation of the FFT algorithm for the group of integers modulo 2" in the quantum 



context, showing how the group-theoretic formalism leads to the standard quantum 
network and identifying the property of entanglement that gives rise to the exponen- 
tial speedup (compared to the classical FFT). Finally we outline the use of the Fourier 
Ch I transform in extracting periodicities, which underlies its utility in the known quantum 

^ ' algorithms. 



Introduction 



In 1982 Feynman[]I]] noted a profound difference in the nature of physical evolution 
5^ \ governed by the laws of quantum physics as compared to evolution under the laws 

of classical physics. He observed that quantum mechanics (apparently) cannot be 
efficiently simulated on a classical computer (or by any classical means) i.e. that the 
simulation of a general quantum evolution by any classical means appears to involve 
an unavoidable exponential slowdown in running time. This observation embodies 
the essence of the subject of quantum computation so we will begin by elaborating 
its meaning in terms of a simple example. 

Consider a discrete sequential quantum process defined as follows. We start with 
a row of qubits (i.e. 2 level systems with a preferred basis labelled {|0),|1)}) all 
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initially in state |0): 



|0) |0) |0) ■■■ |0) ••• 

qubit 1 qubit 2 qubit 3 ■ ■ ■ qubit j • • • 

We are also given a fixed 2-qubit interaction (or 2-qubit "quantum gate" ) U which is 
a unitary operation that may be applied to any selected pair of qubits. Further- 
more we have a program of instructions specifying the pairs of qubits to which the 
gate should be sequentially applied. Thus step k of the program is "apply U to qubits 
{ikiJk) and replace them in the row", for = 1, . . . ,n. After n steps in this process 
we measure qubit 1 in its preferred basis obtaining or 1 according to a probability 
distribution Vn = {Pn(0),Pn(l)}- Thus by implementing this process in an actual 
quantum physical system we can sample the distribution Vn in time 0{n) i.e. after a 
time which grows linearly with n. 

Our problem is to mimick this process by classical means. More precisely, we wish 
to describe a classical probabilistic process which enables us to sample the distribution 
Vn defined by the above quantum process. A simple way of achieving this is the 
following. The operation U is just a 4 x 4 unitary matrix and given the starting 
state with the program, we can sequentially compute by hand - using simple matrix 
multiplication - the quantum state at each successive stage k. Then knowing the 
state at stage n, the rules of quantum measurement theory enable us to calculate 
p„(0) and Pn(l), so finally we toss a correspondingly biassed coin. 

This classical simulation has the following notable characteristic feature. The 
quantum state after k steps is generally a fc-qubit entangled state requiring 0(2'') 
coefficients for its description (as typically an extra qubit may be brought in at each 
step). Thus because of entanglement, we have an exponential growth in time of the 
information needed to describe the state. Hence the classical simulation will slow 
down exponentially in time under the weight of this exponentially growing informa- 
tion that needs to be processed in each step. To sample Vn our classical simulation 
will require 0(2"^) time while the quantum process marches ahead in unflagging lin- 
ear time. There is no more efficient classical method known to solve this problem. 
Thus according to the laws of quantum mechanics. Nature remarkably is able to pro- 
cess information exponentially more efficiently than can be achieved by any classical 
means! 

Note that if only product states of qubits were available then the information 
needed to describe the state would grow only linearly with n (being n times the 
amount of information needed to describe a typical single qubit state). Thus the 
exponential speedup in our example of quantum information processing is fundamen- 
tally a feature of quantum entanglement. This point has been elaborated in Jozsa|p!2|. 
Indeed it provides an extraordinary manifestation of entanglement which is entirely 
independent of the auxiliary notion of non-locality. 
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We may attempt to mimick the quantum process using classical waves which ad- 
mit the possibility of superposition of modes. For example we might represent each 
qubit by a vibrating elastic string with fixed endpoints and select two lowest energy 
modes of vibration to represent the states |0) and |1). It is then possible to construct 
the general superposition corresponding to a |0) + 6 |1). However, regardless of how 
much the strings interact with each other in their subsequent (externally driven) vi- 
brational evolution, their joint state is always a product state of n separate vibrations. 
The total state space of the total classical system is the Cartesian product of the indi- 
vidual state spaces of the subsystems whereas quantum-mechanically, it is the tensor 
product. This crucial distinction between Cartesian and tensor products is precisely 
the phenomenon of quantum entanglement. Nevertheless, we may yet attempt to 
represent entanglement using classical waves in the following manner. The state of n 
qubits is a 2" dimensional space and can be isomorphically viewed as the state space 
of a single particle with 2" levels. Thus we simply interpret certain states of a single 
2" level particle as "entangled" via their correspondence under a chosen isomorphism 
between 0"?^2 and 7^2" (where Tik denotes a Hilbert space of dimension k.) In this 
way, 2" modes of a classical vibrating system can apparently be used mimick general 
entanglements of n qubits. However the physical implementation of this correspon- 
dence appears always to involve an exponential overhead in some physical resource so 
that the isomorphism is not a valid correspondence for considerations of complexity 
i.e. when the amount of physical resources required to achieve the representation 
is taken into account. For example suppose that the 2" levels of the one-particle 
quantum system or corresponding classical system, are equally spaced energy levels. 
A general state of n qubits requires an amount of energy that grows linearly with n 
(since we will need at most to excite each qubit to its upper level) whereas a general 
state of the 2" level quantum or classical system requires an amount of energy that 
grows exponentially with n. To physically realise a system in a general superposition 
of 2^ modes we need exponential resources classically and linear resources quantum 
mechanically because of the existence of entanglement. 

Our discussion above about the information needed to describe a state, indicates 
that n qubits have an exponentially larger capacity to represent information than 
n classical bits. Note that although n classical bits have 2"' possible states, each of 
these states may be described by just n bits, in contrast to the quantum situation 
where 0(2") superposition components may be involved in a single state. However the 
information embodied in the quantum state has a further remarkable feature - most of 
it is inaccessible to being read by any possible means! Indeed quantum measurement 
theory places severe restrictions on the amount of information that we can obtain 
about the identity of a given unknown quantum state. This intrinsic inaccessibility of 
the information may be quantified P, Q in terms of Shannon's information theory ^j. 
In the case of a general state of n qubits, with its 0(2") information content, it turns 
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out that at most n classical bits of information about its identity may be extracted 
from a single copy of the state by any physical means whatsoever. This coincides 
with the maximum information capacity of n classical bits. 

The full (largely inaccessible) information content of a given unknown quantum 
state is called quantum information. Natural quantum physical evolution may be 
thought of as the processing of quantum information. Thus the viewpoint of compu- 
tational complexity reveals a new bizarre distinction between classical and quantum 
physics: to perform natural quantum physical evolution, Nature must process vast 
amounts of information at a rate that cannot be matched by any classical means, yet 
at the same time, most of this processed information is kept hidden from us! However 
it is important to point out that the inherent inaccessibility of quantum information 
does not cancel out the possibility of exploiting this massive information processing 
capability for useful computational purposes. Indeed, small amounts of information 
may be extracted about the overall identity of the final state which would still require 
an exponential effort to obtain by classical means. The ability to sample the prob- 
ability distribution Vn above provides an example. A more computationally useful 
example is given by the technique of "computation by quantum parallelism" ||^, ^ 
according to which a superposition X]?=i of exponentially many input values i for 
a function / may be set up in linear time and a single subsequent function evaluation 
will provide exponentially many function values in superposition as |/) = X^i N) |/(^))- 
The full quantum information of this state incorporates the information of all the in- 
dividual function values f{i) but this is not accessible to any measurement. However 
certain global properties of the collection of all the function values may be deter- 
mined by suitable measurements on |/) which are not diagonal in the standard basis 
{\i) \j)}. For example if / is a periodic function, we may determine the value of the 
period [0, which falls far short of characterising the individual function values but 
would generally still require an exponential number of function evaluations to obtain 
reliably by classical means. 

Entanglement Enhanced Information Processing 

Suppose that we have a physical system of n qubits in some entangled state {ip) and 
we apply a 1-qubit operation U to the first qubit. This would count as one step 
in a quantum computation (or rather a constant number of steps independent of n, 
if U needs to be fabricated from other basic operations provided by the computer). 
Consider now the corresponding classical computation, {tp) may be described in com- 
ponents (relative to the product basis of the n qubits) by aij...j„ where each subscript 
is or 1, and U is represented by a 2 x 2 unitary matrix U- . The application of U 
corresponds to the matrix multiplication 
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Thus the 2x2 matrix muhiphcation needs to be performed 2"^^ times, once for 
each possible value of the string 22 ■ ■ ■in, requiring a computing effort which grows 
exponentially with n. On a quantum computer, because of entanglement, this 2"~^ 
repetition is unnecessary. 

Consider now a unitary transformation U of n qubits (or more precisely a family 
of such transformations labelled by n). U may be described by a 2" x 2" matrix and 
the computation of U classically by direct matrix multiplication requires 0(2"2'^) 
operations. Even on a quantum computer U needs to be fabricated ( "programmed" ) 
out of the basic operations provided by the computer, each of which operate only on 
some constant number of qubits. In general U will require an exponential number of 
such basic operations for its implementation. It may be shown || that 0(2"2") 
operations will always suffice to program U to any desired accuracy. 

Suppose now that U has the following special form. Let c be any constant, inde- 
pendent of n. Suppose that U consists of the sequential application of p{n) unitary 
operations Vi,i = 1, . . . ,p{n) where each Vi operates on only some c out of the n 
qubits and p{n) is a polynomial in n. An immediate generalisation of the argument 
above shows that each Vi may be classically implemented (by matrix multiplication) 
in 0(c^2"~'^) = 0(2'*) steps so that the classical computation of U now requires 
0(p(n)2") steps. This represents an exponential saving over a general U which re- 
quired 0{2^2^) steps but it is still exponential in n. An important example of this 
partial exponential speedup for classical computation is the so-called fast Fourier 
transform algorithm [1^, as compared to the regular Fourier transform algorithm. 



On a quantum computer each Vi requires some constant (independent of n) number 
of steps to implement (programming the c-qubit operation Vi in terms of the basic 
operations) so that U requires only p{n) steps to implement. In summary, if U has the 
special form given above then it still requires exponential time to compute classically 
(although it does provide a partial exponential benefit here already) but it requires 
only polynomial time to compute on a quantum computer. Note however that after 
the quantum computation only a small amount of information about the transformed 
data is accessible to measurement, whereas the classical computation allows the full 
information to be accessed. 

The Super-fast Quantum Fourier Transform 

The Fourier transform on a finite Abelian group G is a large unitary operation which 
arises naturally in the mathematical formalism of group representation theory. Fur- 
thermore it factorises in the special way described in the previous section if the 
group has some additional structure and it is known to be a basic tool for vari- 
ous useful computational tasks, in particular the problem of determining periodicity. 
Consequently, in view of the discussion above, it can lead to quantum algorithms 
|T3| , H, H, H, |10|, |Tl], [iSl which run substantially faster than any known classical 
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algorithm for the corresponding computational task. In this section we will outline 
the construction of the Fourier transform and describe its factorisation into unitary 
operations of a constant size. 

Let {G, +) be any finite Abelian group where we write the group operation in 
additive notation. Let |G| denote the number of elements of G. An irreducible 
representation of G is a function 

X-.G^C* 

(where C* denotes the non-zero complex numbers) satisfying 

X{9i + 92) = X{9i)x{92) (2) 

i.e. X is a group homomorphism from the additive group G to the multiplicative 
group G*. The condition eq. (0) has the following consequences (see e.g. |2^, [13| for 
proofs) . 

(A) Any value x(fl') is a \G\^^ root of unity. Thus x i^ay be viewed as a group 
homomorphism x '■ G ^ where is the circle group of all unit modulus 
complex numbers. 

(B) Orthogonality (Schur's lemma): If Xi and Xj are any two such functions then: 

Yl Xi{9)Xj{9) = Sij (3) 
l<^l gee 

(where the overline denotes complex conjugation). 

(C) There are always exactly |G| different functions x satisfying eq. (|). 

In view of (C) these functions may be exhaustively labelled by the elements of G. 
Let {xg '■ 9 G G} be any such chosen labelling. Then the Fourier transform on G is 
the |G| X \G\ matrix JF whose rows are formed by listing the values of the functions 



y^gk = ^=Xg{k) g,keG (4) 
\G\ 



Xg- 



Note that by (B) JF is always a unitary matrix. 

In the context of quantum computation we will have a Hilbert space Ti. of dimen- 
sion \G\ with a basis {\g) : g E G} labelled by the elements of G. Thus there is a 
natural shifting action of G on Ti given by 

U{k):\g)-^\g + k) k,geG (5) 
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These operations all commute since G is Abelian so there exists a basis of simultaneous 
eigenstates of all the shifting operators. According to (B) the states 



\Xk) = ^ E Xk{g) \g) keG (6) 

\ \G\ geG 



form an orthonormal basis of Ti and using eq. (H) we get 

U{g) \xk) = e^'=(^) \xk) 

so that {\Xg) : (? G G} is the basis of common eigenstates of the shift operators. This 
basis is also called the Fourier basis. The Fourier transform JF is a unitary operation 
on H and using eq. @ with (H) and property (B) we readily get: 

^\Xs) = \9) (7) 
so that the Fourier transform interchanges the standard and Fourier bases. 



Let Zq denote the additive group of integers mod q. It is well known [|1J] that any 
finite Abelian group G is isomorphic to a direct product of the form 

G = X X ... X Zmr (8) 

(Furthermore we may require that rrii divides mj+i and then the numbers rrii are 
unique). If we assume (usually without loss of generality) that the group G is pre- 
sented as a product of the form eq. (§), then we can explicitly describe the irreducible 
representations (0) and obtain a canonical labelling of them by the elements of G. 
Suppose first that G = Zm- Consider the group homomorphism given by 



r: GxG S 
(a, 6) — i> e 



It is easily verified that for each fixed a E G the function Xa '■ G ^ given by 
Xa{b) = T{a,b) satisfies eq. @ and there are \G\ such functions. Thus we have 
obtained an explicit formula for the irreducible representations, labelled in a natural 
way by the elements of G. For the general case of a product G = Z^x x Zm^ x . . . x Zm^ 
we simply multiply the corresponding factors in eq. (||) obtaining 



r : G X G ^ 

((ai,...,a,),(6i,...,6.)) ^ exp27r^(^ + ^^--- + ^ 



(10) 



rrtr 



and again 

Xg,{9i) = r{gx,g2) (11) 
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provides the irreducible representations labelled by the elements of G. 

As an example consider the group (^'2)" of all n-bit strings. From eqs. (p!OD (p 
and (I) we see that the Fourier transform is just 



V2" 



where a-p = siti + . . . + s„t„ mod 2 if cr = Si . . . s„ and z/ = ti . . . t„. Thus in this case 
the Fourier transform coincides with the Hadamard (Walsh) transform. If G = 



then we see using eqs. (H) and (|TTD that 

1 



= -—e^^'^ a, 6 = 0, ... 2" 



giving the familiar discrete Fourier transform modulo 2". 

As a unitary matrix the Fourier transform will act on vectors of length \G\. We 
may view any such vector as a function f : G ^ C on G whose list of values 
f{gi), . . . , f{g\G\) defines the vector. The Fourier transform of / is then given by 

m = E ^kjig) = ^ E Xk{g)f{g) k e G (12) 



We now describe the basic factorisation property of this large unitary transforma- 
tion which is necessary for its efficient (i.e. polynomial time) implementation in the 
context of quantum computation. The factorisation will be carried out relative to a 
subgroup H of G and again the key ingredient will be the property given by eq. (|). 
The basic technique was developed by Cooley and Tukey leading to the so-called 



fast Fourier transform (FFT) algorithm in classical computation (which provides the 
partial exponential speedup noted in the previous section) but the essential idea oc- 
curs already in the work of Gauss |p!6 |. 



Let if be a subgroup of G with index / = IGl/lifl. Let ki+H, k2+H, . . . , kj+H be 
a complete list of the cosets of where k+H ^ G denotes the subset given by {k+h : 
h G H}. Thus G is partitioned as a disjoint union (fci + if) U (/c2 + //) U . . . U (fc/ + H). 
Hence the elements g & G may be written in a unique way in terms of the cosets as 
g = ki + h. Using eqs. (|12D and (|^) we get: 



1 _ „, , . . I ' 



f{i) = J2fi9)xi{g) = ^j:J2fik^ + h)xi{h + h) 



1 



,ExKfci)E/^Wx/W (13) 

G\ i=i heH 
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where fi for i = 1, . . . , I are the functions on H defined by the restrictions of / to the 
cosets: fi{h) = f{ki + h). The functions xi restricted to the subgroup H satisfy eq. 
@ on H so they are irreducible representations of H . Hence the sum over H in eq. 
( p!3D amounts to evaluating the Fourier transform on H of the functions /j. Thus eq. 
( p!3D expresses a decomposition of the Fourier transform on G into the evaluation of I 
Fourier transforms on H whose results are then combined linearly in sums of length I 
with coefficients Xii^i): done for each / e G. Hence the number of operations required 
is 

Oi\H\^ xI+\G\xI) = 0{\G\i\H\ + I)) (14) 

where we have used / = IGl/l-ffl. This is generally better than the 0(|G|.|G|) oper- 
ations for the direct (matrix multiplication) calculation of the Fourier transform on 
G. For example we choose H so that / is small, say 1 = 2 giving \H\ = \G\/2 and 
then eq. (|1^ represents an approximate halving of running time. 

To enhance this benefit we iterate the construction on a tower of subgroups 

Gd HiD H2D ■■■D HnD {0} 

of greatest possible length, ultimately expressing the Fourier transform of G in terms 
of that on the (small) subgroup Hn- An extensive survey of this technique is given 
in 113. We will illustrate it here only for the g and discuss the effect of the 

resulting decomposition on the quantum computational implementation. has an 
optimal tower of subgroups with each successive inclusion having the minimal possible 
index of 2: 

Z^"- D Z2n-\ D Z2n,-'2. D ■ ■ • Z> Z2 3 {0} 

(Here ^2^-1 is the subgroup {0,2,4, ... ,2" — 2} of all even integers in ^2", Z2n-2 is 
the subgroup {0, 4, 8, . . .} of all multiples of 4 etc. and Z2 is the subgroup {0, 2"^^} 
). Consider a general position Z^^ D ^2^-1 in this chain and let FT2™ denote the 
Fourier transform Oil ZjOm . The irreducible representations of Z2 

Xj{k) = {w^f for j,fc = 0,...2'"-l (15) 

where w = exp Then eq. ( [TB| ) becomes (writing out the i-sum explicitly): 



k=0 V ^ 



1 2jk 2—^-1 2jk 



(16) 



Here the f{2k) in the first sum and f{2k + 1) in the second sum give the function / 
restricted respectively to the cosets of Z2m-i C ^2™ (i.e. the even and odd positions 
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Ill ZjOm, ). Note that the irreducible representations of Z2 

m--l cirG the functions eiven 

in eq. (|T5|) with w replaced by vP' = exp^^. Thus the two fc-sums on RHS of 
eq. ( |T6|) are just FT2m-i of the even and odd labelled values of /. As j in eq. (|16D 
runs through the values to 2™ — 1, we cycle twice through the 2™~^ components 
of the FT2m-i's (noting that (w^)^™ = !)• If we restrict j to running through 
the values to 2™"^ — 1 then /(j) and f{j + 2*^"^) are both obtained from the j^^ 
components of the two FT2m-i transforms on RHS of eq. ([T6|) , combined respectively 
with coefficients ■^{l,w^) and ■^{l,w^~^^"' ) = -^{1, —w^). Thus eq. (|TB|) may be 
described as 

/(j) = 75( J*' Cpt. of FT(/e.en) + ■ J*^ Cpt. of FT(/odd) ) 1 

f{j + 2—1) = ^( cpt. of FT(/even) " ■ Cpt. of FT(/odd) ) J 

where /even and /odd refer respectively to the 2™'"^ even and odd labelled values of / 
and j ranges from to 2"*"^ — 1. Now if C(2™') denotes the number of operations 
required to (classically) compute FT2m then eq. ([16D shows that 

^(2™) = 2C(2™-i) + 0(2™) 

where the 0(2™") arises from the extra additions and multiplications needed for the 
2"^ j-values in eq. ([iBD , to linearly combine the results of the two FT2m-i operations. 
The solution of this recursion relation is 

C(2") = 0(722'^) 

giving the partial exponential speedup (compared to 0(2^^2")) noted previously. 

In the context of quantum computation the data values /(j) for j = 0, . . . , 2™ — 1 
reside in the amplitudes of an entangled state |/) of m qubits. Writing j in binary 
as an m bit string we have 

1 

I/) = l/(im-l---jljo)) bm-l)---|jl) lio) 

jo JlvJm-l=0 

and the qubits are numbered 0, 1, . . . , m — 1 from right to left. The two FT2m-i opera- 
tions in eq. (|16|), which operate on even and odd numbered components respectively, 
may then be implemented by a single FT2m^i operation on qubits m — l,m — 2, 
since the values and 1 of the remaining rightmost index respectively determine the 
even and odd labelled positions (c.f. the discussion of eq.(|l|)). The j^^ component of 
FT2m-i (/even) (respectively FT2m-i(/odd)) then resides as the amplitude in dimension 
2j (respectively 2j + l). Thus to perform the linear recombination of the two FT2m-i's 
eq. (0) shows that we need to 
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(a) perform the unitary operation 



V2 [ 1 -w^ ) 

on dimensions (2j, 2j + 1) for each j = 0, . . . , 2™"^ — 1. 

(b) Reorder the answers according to the permutation (2j, 2j + 1) — > + 2™"^) 
for each j = 0, . . . , 2™"^ — 1 to get /(j) as the amphtude in dimension j. 

This would appear to involve exponentially many operations (for the 2™"^ values of 
j) but using the entanglement effects discussed at eq. (||), we can achieve the result 
with only 0{m) operations as follows. Note first that 

^[l -w^ ) ^[l -1 )[0 J 

The operation Bj in dimensions (2j, 2j ' + 1) leaves the even dimension unchanged 
and applies an phase shift in the odd dimension. This may be achieved for all 
j values simultaneously by applying a 2-qubit gate Cp to qubits 1 and p for each 
p = l,...,m — 1. Here Cp is the conditional phase shift of w^'' applied to qubit p 
only if both qubits and p are 1. In the standard basis of qubits and p we have: 



/I 











\ 





1 
















1 







V 








op-l 


) 



Using the entanglement effects described at eq. (|l|) we see that the successive appli- 
cation of the m — 1 operations Cp builds up a phase of in dimension 2j + 1 for 
each j. The requirement that the zeroth qubit have value 1 selects the odd positions 
and the conditional phase shift in Cp builds up the value successively for each '1' 
in the binary expansion of j. All (exponentially many) values of j with '1' in the 
p^^ place are treated simultaneously. Finally the 1-qubit operation H is applied just 
once to qubit 0, which simultaneously applies H to all pairs (2j, 2j ' + 1) given by all 
possible values of the remaining indices for qubits 1 to m — 1 (c.f. eq. (|l|)). 
To implement (b) i.e. the permutation of dimension labels given by 

even labels: 2j j odd labels: (2j + 1) ^ (j + 2""-^) 

we simply cyclically permute the qubit labels as m-bit strings: 

im-i ■ ■ ■ hio — ^ ioim-i ■ ■ - ii 
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If the label was even (i.e. io = 0) then the value is halved and if it was odd {io = 1) 
then the cycling of io to the leading position adds 2"^"^ and the residual even part is 
halved. This cycling may be physically achieved by m — 1 state swaps, of qubits 
and 1, then 1 and 2 etc. up to qubits m — 2 and m — 1. Alternatively we may just 
reorder the output wires as shown in figure 1 below. 








EQUALS: 



■ m 




Figure 1. The network diagram for the decomposition of FT2m into FT2m-i and 
0(m) extra operations. The conditional phase shift Cp on the p*^ qubit is denoted by 
a box on the p^^ qubit line and a connection across to the O*'^ qubit with a cross to 
denote the fact that its operation on the p^^ qubit is "controlled" by the requirement 
that the 0*"^ qubit have value 1. 



Iterating this construction for FT2m-i in terms of FT2m-2 etc. yields the standard 
network for the fast Fourier transform on as given for example in ||10|| . 

If (5(2™") denotes the number of operations needed to implement FT(2™') in the 
quantum context then the above description shows that 

Q{2^)=Q{2^-')+0{m) 
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giving 

Q{T) = 0(n2) 

This quadratic time quantum algorithm for FT(2") is used in Shor's factoring algo- 
rithm rui, 111. 



Utility of the Fourier Transform 

The utility of the Fourier transform JF in the algorithms of Deutsch, Simon and Shor 



Ig, I], ^ has been described in ||T3[. We will here outline in general terms, its funda- 
mental application to the determination of periodicities. A different interpretation of 
T in terms of the problem of phase estimation, has been given in 

Let / : G ^ X be a function on the group (taking values in some set X) and 
consider 

K = {keG: f{k + g) = fig) for all g G G} 

K is necessarily a subgroup of G called the stabiliser or symmetry group of /. It 
characterises the periodicity of / with respect to the group operation of G. Given 
a device that computes /, our aim is to determine K. More precisely we wish to 
determine K in time 0(poly(log |G|)) where the evaluation of / on an input counts as 
one computational step. (Note that we may easily determine K in time 0(poly(|G|)) 
by simply evaluating and examining all the values of /). We begin by constructing 
the state 

l/) = ^El^?)l/(^?)) 

\J\G\ geG 

and read the second register. Assuming that / is suitably non-degenerate - in the 
sense that f{gi) = f{g2) gi — g2 & K i.e. that / is one-to-one within each period 
- we will obtain in the first register 

mgo)) = E \9o + k) (18) 



fcex 

corresponding to seeing /((^o) in the second register and g^ has been chosen at random. 
In eq. ([T8|) we have an equal superposition of labels corresponding to a randomly 
chosen coset of K in G. Now G is the disjoint union of all the cosets so that if we 
read the label in eq. ([I8| ) we will see a random element of a random coset, i.e. a 
label chosen equiprobably from all of G, yielding no information at all about K. The 
Fourier transform will provide a way of eliminating g^ from the labels which may 
then provide direct information about K. Consider the basis {\Xg) : (? G G} of shift 
invariant states introduced in eq. (|^). Next note that the state in eq. (^8|) may be 
written as a ^fQ-shifted state: 

E \9o + k) = 
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Hence if we write this state in the basis {\Xg) i9 ^ G} then 1^) and X^fc Is'o + k) will 
contain the same pattern of labels, determined by the subgroup K only. According 
to eq. (0) the Fourier transform converts the shift-invariant basis into the standard 
basis. Thus after applying to eq. (p!8| ) we may read the shift-invariant basis label 
by reading in the standard basis, yielding information about K. 

In terms of the presentation of G given in eq. (^) and the associated formulas for 
the irreducible representations given by eqs. (0) and ( |TT|) we may compute explicitly 
the pattern of labels associated with a subgroup K G G. As an example consider 
G = Zmn and K = mZ = {0, m, 2m, . . . ,{n — l)m} with \K\ = n. Then the Fourier 
transform of the fundamental periodic state \K) = -^J2keK \k) is 

^\^) = Z^T.\T.Xi{k)] 10 (19) 



Thus the labels appearing are precisely those / G Zmn for which 

E Xiik) ^ (20) 

k(^K 

To sort out this condition we introduce a further elementary property of irreducible 
representations. For any group G the constant function xio) = 1 fo^' all g G G, 
is clearly an irreducible representation (the trivial representation) and using the or- 
thogonality property (B) between x and any other irreducible representation x' we 
see that 

E x'ig) = 

Now xi restricted to the subgroup K is an irreducible representation of K so eq. (pO|) 
can hold if and only if 

Xi{k) = 1 for all k e K 

According to eqs. ([ToD and (0) we have 

/ 7 \ kl ^ -cl 

Xi\k) = exp27r« = exp27r2— 

mn n 

where we have introduced c using the fact that k = cm is always a multiple of m, 
by definition of K. This will equal 1 for all c = 0, . . . , (m — 1) if and only if / is a 
multiple of n i.e. Z = 0, n, 2n, . . . , (m — l)n. Thus the pattern of labels associated 



with mZ C Zmn is nZ and furthermore in eq. ([T9| ) each such label will appear with 
equal amplitude A similar calculation for the subgroup {0,^} C (^2)" (where ^ 
is a chosen n-bit string) shows that the resulting pattern of labels, after applying the 
Fourier transform for (.^2)" to the periodic state :^(|0) + |^)), is {u : ■ u = 0}. This 
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fact forms the basis of Simon's algorithm 



PI 



Conclusion 

Let \ip) be an n-qubit entangled state and U a 1-qubit unitary operation. We have 
seen that the one-step physical operation of applying U to (say) the first qubit of 

corresponds to a state transformation which generally requires an exponential (in 
n) effort to compute classically. Indeed mathematically the transformation is repre- 
sented by a tensor product U ® I2® ■ ■ ■ ® h (where I2 is the 2 by 2 identity matrix, 
which represents the operation of "doing nothing" on the corresponding qubits). The 
tensor product spreads the effect of U into an exponentially large matrix. Stated 
otherwise, we can say that the physical operation of doing nothing to a subsystem of 
an entangled system is a highly nontrivial operation and gives rise to an exponen- 
tially enhanced information processing capability (when performed in conjunction 
with some operation on another small part of the system). 

We have given an analysis of the implementation of the fast Fourier transform al- 
gorithm in a quantum context and shown that its exponential speedup (as compared 
to the corresponding classical computation) derives wholly from the above tensor 
product property. We have also given a general discussion of the role of entanglement 
in quantum computation and the utility of the Fourier transform in the known quan- 
tum algorithms. 
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